Huygens&#39; box methodology for signal integrity analysis

ABSTRACT

A method for performing a signal integrity analysis on an integrated circuit (IC) that includes a plurality of scatterers by dividing the scatterers into subgroups using a nested Huygens&#39; equivalence principle algorithm and solving a set of equations realized thereby with a reduced coupling matrix. The method includes decomposing the IC design into a plurality of small non-overlapping circuit sub-domains, wherein each of the sub-domains is formed as a small, enclosed region. Each sub-domain is analyzed independently of the other sub-domains using only electric fields to represent the interactions of each sub-domains with the other sub-domains as equivalent currents on equivalent surfaces of the plurality of sub-domains. Neighboring equivalent sub-domains are grouped together to form larger sub-domains using equivalent currents on equivalent surfaces to represent the interactions of the sub-domains. The steps of analyzing and grouping the sub-domains are repeated until the grouping approaches a box comprising the entire domain, and that the domain interactions between every sub-domain have been analyzed.

BACKGROUND OF THE INVENTION

The invention relates to general electromagnetic simulation and electromagnetic field data processing for signal integrity analysis. More particularly, the invention relates to a method and computer system for solving complicated electromagnetic simulation problems in high-speed integrated circuit (IC) signal integrity analysis.

During electrical analyses of devices that are small with respect to wavelength, such as electronic packaging structures, chips, interconnects, and printed circuit boards, linear systems of equations must be solved. The present invention enhances existing electromagnetic simulation EDA tools using a novel Huygens' box methodology. The novel method and system requires use of only electric field integral equations (EFIE) instead of both EFIE and magnetic field integral equations (MFIE).

A Huygens' box approach to electromagnetic analysis of electronic structures includes constructing a large structure or domain from a set of smaller, non-overlapping sub-domains. The field property of the internal sub-domains are replaced with equivalent sources on the outer surfaces, i.e., discretized point sources. Each sub-domain is analyzed separately, without regard to the others. The method may be used enable conventional EDA tools to analyze complicated structures by compartmentalizing the problem into smaller equivalent units with fewer unknowns. As such, the present invention may serve as an industry standard to organize outputs from different electromagnetic vendor tools, generate parameterized model library databases, and in a readily parallelizable process.

A nested equivalence principle algorithm (NEPAL) that divides scattering objects into subgroups is known. W. C. Chew and C. C. Lu, “NEPAL—An Algorithm for Solving the Volume Integral Equation,” Microwave and Optical Technology Letters, vol. 6, no. 3, pp 185-188, Mar. 5, 1993 (“Chew and Lu (1)”). For each subgroup, the NEPAL algorithm substitutes interior scatterers with boundary sub-scatterers using Huygens' equivalence principle. The subgroups are combined by levels to achieve a computational complexity of O(N²) for non-iterative solutions to 3D scattering problems. The NEPAL Algorithm, however, uses spherical point sources as the equivalent sources so that solving electromagnetic scattering of 3D objects using spherical wave manipulation is very complicated. W. C. Chew and C. C. Lu, “The Use of Huygens' Equivalence Principle for Solving the Volume Integral Equation of Scattering,” IEEE Trans. Antennas and Propagat., vol. 41, no. 7, pp. 897-904, July 1993 (“Chew and Lu (2)”); and C. C. Lu and W. C. Chew, “The Use of Huygens' Equivalence Principle for Solving 3-D Volume Integral Equation of Scattering,” IEEE Trans. Antennas and Propagat., vol. 43, no. 5, pp. 500-507, May 1995 (“Lu and Chew (3)”).

Known alternative techniques include the use of a combination of surface integral equations, and apply a form of Huygens' equivalence principle. M.K. Li and W.C. Chew, “Wave-Field Interaction with Complex Structures Using Equivalence Principle Algorithm,” UIUC CCEM Lab Research Report, 2006 (“Li and Chew (4)”), and M. K. Li, W. C. Chew and L.J. Jiang, “A Domain Decomposition Scheme to Solve Integral Equations Using Equivalent Surfaces,” 2006 IEEE Antennas and Propagation Society International Symposium, pp. 2897-2900, Albuquerque, N.M., USA (“Li and Chew (5).”). These alternative techniques use both electric field and magnetic field to build the field response for an enclosed sub region to mimic the general [S] parameters. Boundary equivalent sources based on surface basis are used to preserve the sub region interactions. A nesting process similar to NEPAL is employed to combine sub regions. Since both EFIE (Electric Field Integral Equation) and MFIE (Magnetic Field Integral Equation) are involved, however, many operators are needed, rendering such an implementation process complicated.

Also known is a reduced-coupling method implemented to deduce the dense interactions through boundary ports. B. J. Rubin, A. Mechentel, and J. H. Magerlein, “Electrical Modeling of Extremely Large Packages,” 1997 Electronic Components and Technology Conference, pp 804-809 (“Rubin (6)”); and US Pending Application Serial No. 2005/0222825 to Feldmann, et al., Method, Apparatus and Computer Program Providing Broadband Preconditioning Based on Reduced Coupling For Numerical Solvers (“the '825 application”). Rubin (6) and the '825 application accelerate the computational process, but not rigorously since the connecting boundaries are not closed.

OBJECTS AND SUMMARY OF THE INVENTION

In view of the shortcomings of the prior art, it is an object of the invention to realize a systematic mathematical methodology to partition electromagnetic simulation into smaller jobs that can be run much faster, using a less complicated technique or method.

It is another object to realize a methodology to join different electromagnetic field simulation processes, which processes maintain the physically correct wave-properties of any structure without any simplification.

It is another object of the invention to realize a physically and mathematically rigorous approach to replace conventional [S] parameter approximation currently widely used in the industry to solve scattering problems, and extend its application to much higher operating frequencies based on wave properties.

The invention is anticipated to be applied as an industry standard for pre-physical design (pre-PD) stage to unify outputs from different electromagnetic tools. The inventive method conveniently improves electric field integral equation (EFIE) based solvers for operation upon much more complicated structure simulations, providing a powerful tool to build customized libraries of models for re-use and for parameterization of general chip interconnects, and packaging design and analysis, as well as a simulation framework for fast computational algorithms including FFT, Fast Multipole Algorithm, and multi-grid method, etc. The invention further provides a methodology for use as a fast-preconditioner for iterative solutions, and readily parallelized on multiple server systems.

DESCRIPTION OF THE DRAWING FIGURES

FIG. 1 is a depiction of direct coupling interaction between two perfect electrical conductor (PEC) objects;

FIG. 2 is a depiction of the PEC objects of FIG. 1 wherein PEC I is enclosed in a Huygens' Equivalence Box, where PEC II is not; and

FIG. 3 is a flowchart highlighting a whole solution domain decomposed into subdomains in accordance with a novel equivalence principle algorithm (EPA) of the invention.

DETAILED DESCRIPTION OF THE INVENTION

The invention includes a novel, simple and convenient methodology to implement the Huygens' Equivalence Principle algorithm using only electrical field integral equations (EFIE). The novel method may be treated as an extension of the reduced coupling method for implementing a rigorous solution. The novel method implements a variation on the NEPAL nesting scheme described in detail in Chew and Lu (1), Chew and Lu (2) and Lu and Chew (3), the Equivalence Principle Algorithm disclosed and described in Li and Chew (4) and Li and Chew (5), and implementation of reduced coupling method as described Rubin (6) and the '825 application.

The Huygens' Box Approach can be thought of as a mechanism for constructing a large structure, intended for electromagnetic analysis, from a set of smaller, non-overlapping regions. The Huygens' box approach is analogous to building a large model from Lego blocks, wherein each block is analyzed separately without regard to the others. That is, the equivalence boxes separate and enclose analyzed structures into small closed regions. The equivalent electric current and magnetic current on the EPA surfaces are formulated together with the realistic electric current and magnetic current on the conductor and dielectric surfaces. The blocks are combined through a mathematical procedure to reform the original structure. The interactions between different regions are forced to go through the closed surfaces of the boxes. A direct mapping of this interaction decomposition to a numerical model is defined by a set of surface integral equations, which are EFIE only in this method.

FIGS. 1 and 2 together highlight the coupling interaction between perfect electrical conductors (PECs) based on the Huygens' equivalence principle. FIG. 1 depicts direct coupling interaction or the scattering between two PEC objects, where FIG. 2 depicts the first PEC (PEC I of FIG. 1) enclosed in a Huygens' equivalent box. With PEC I enclosed as shown in FIG. 2, and where PEC II is not enclosed in a Huygens equivalent box, the resultant MOM matrix is arranged as follows:

$\begin{matrix} {{\begin{bmatrix} L_{11} & 0 & L_{11^{\prime}} & K_{11^{\prime}} \\ 0 & L_{22} & L_{21^{\prime}} & K_{21^{\prime}} \\ L_{1^{\prime}1} & 0 & L_{11^{\prime}} & K_{11^{\prime}} \\ 0 & L_{1^{\prime}2} & L_{11^{\prime}} & K_{11^{\prime}} \end{bmatrix}\begin{bmatrix} J_{1} \\ J_{2} \\ J_{1^{\prime}} \\ M_{1^{\prime}} \end{bmatrix}} = \begin{bmatrix} 0 \\ {- E_{2}^{inc}} \\ {- E_{1^{\prime}}^{inc}} \\ {- E_{1^{\prime}}^{inc}} \end{bmatrix}} & {{Matrix}\mspace{14mu} (1)} \end{matrix}$

Since the matrix is highly sparse, numerous sparse matrix techniques can be used to optimize the process for solving Matrix (1). W. H. Press, S. A Teukolsky, W. T. Vetterling, and B. P. Flannery, “Numerical Recipes in C++: The Art of Scientific Computing,” 2^(nd) Edition, Cambridge University Press, 2002. Because each of the coefficients has explicit physical meanings, Matrix (1) can be further reduced to a matrix problem in which the details inside the Huygens' box may be eliminated. Using the compression process, a much smaller matrix equation may be arranged as:

$\begin{matrix} {\begin{bmatrix} L_{22} & L_{21^{\prime}} & K_{21^{\prime}} \\ 0 & {L_{1^{\prime}1^{\prime}} - {L_{1^{\prime}1}L_{11}^{- 1}L_{11^{\prime}}}} & {K_{1^{\prime}1^{\prime}} - {L_{1^{\prime}1}L_{11}^{- 1}K_{11^{\prime}}}} \\ L_{1^{\prime}2} & L_{1^{\prime}1^{\prime}} & K_{1^{\prime}1^{\prime}} \end{bmatrix}{\quad{\begin{bmatrix} J_{2} \\ J_{1^{\prime}} \\ M_{1^{\prime}} \end{bmatrix} = \begin{bmatrix} {- E_{2}^{inc}} \\ {- E_{1^{\prime}}^{inc}} \\ {- E_{1^{\prime}}^{inc}} \end{bmatrix}}}} & {{Matrix}\; (2)} \end{matrix}$

If both conductors PEC I and PEC II are enclosed inside different Huygens boxes, the matrix solution becomes:

$\begin{matrix} {\begin{bmatrix} L_{1^{\prime}1^{\prime}} & K_{1^{\prime}1^{\prime}} & L_{1^{\prime}2^{\prime}} & K_{2^{\prime}2^{\prime}} \\ {L_{1^{\prime}1^{\prime}} - {L_{1^{\prime}1}L_{11}^{- 1}L_{11^{\prime}}}} & {K_{1^{\prime}1^{\prime}} - {L_{1^{\prime}1}L_{11}^{- 1}K_{11^{\prime}}}} & 0 & 0 \\ L_{1^{\prime}1^{\prime}} & K_{1^{\prime}1^{\prime}} & L_{2^{\prime}2^{\prime}} & K_{2^{\prime}2^{\prime}} \\ 0 & 0 & {L_{2^{\prime}2^{\prime}} - {L_{2^{\prime}2}L_{22}^{- 1}L_{22^{\prime}}}} & {K_{2^{\prime}2^{\prime}} - {L_{2^{\prime}2}L_{22}^{- 1}K_{22^{\prime}}}} \end{bmatrix}{\quad{\begin{bmatrix} J_{1^{\prime}} \\ M_{1^{\prime}} \\ J_{1^{\prime}} \\ M_{2^{\prime}} \end{bmatrix} = \begin{bmatrix} {- E_{1^{\prime}}^{inc}} \\ {- E_{1^{\prime}}^{inc}} \\ {- E_{2^{\prime}}^{inc}} \\ {- E_{2^{\prime}}^{inc}} \end{bmatrix}}}} & {{Matrix}\mspace{20mu} (3)} \end{matrix}$

which can be directly calculated based on Matrix (2). Therefore, if we want to discard the unknowns inside a Huygens' box, we can modify the interaction between the objects and the box using the following matrix equation coefficients:

$\begin{matrix} {\left. L_{i^{\prime}i^{\prime}}\Rightarrow{L_{i^{\prime}i^{\prime}} - {L_{ii}L_{ii}^{- 1}L_{{ii}^{\prime}}}} \right.\left. K_{i^{\prime}i^{\prime}}\Rightarrow{K_{i^{\prime}i^{\prime}} - {L_{ii}L_{ii}^{- 1}K_{{ii}^{\prime}}}} \right.} & {{Matrix}\mspace{14mu} (4)} \end{matrix}$

Matrix (4) is derived by discarding the columns and rows related to the do-not-care unknowns. Matrix (2) and Matrix (3) are related to Lu and Chew (3) and Li and Chew (5), but obtained using different processes and therefore different operators.

The inventive method requires only using electric field instead of both electric and magnetic field calculations, resulting in a solution and recombination of the separate solutions in the reduced coupling matrix form as described in the pending '825 application. As a result, many interactions from the region boundaries are eliminated compared to the conventional direct methods. By eliminating region boundaries, formulation of the same physical process becomes much simpler and therefore faster. The corresponding storage becomes smaller, and solution speeds increase. The Huygens' Box nesting provides the overall strategy to solve complicated structures hierarchically, preserving the true field properties of the structure, unlike conventional [S] parameter based subdivisions. S parameters only interact through a few ports while the entire surface of the Huygens' Box preserves the field properties and interactions to the adjacent box.

In FIG. 3, the whole solution domain, such as a large interconnect network, electrical packages, etc., is shown decomposed into several sub-domains. Each sub-domain is solved independently at the lowest level and the solution is represented as equivalent currents on equivalent surfaces. Then, the neighboring equivalent surfaces are grouped together to form a larger box. The representation of the electromagnetic properties of the sub-domains includes equivalent currents that only reside on the equivalent surfaces. The process moves forward iteratively until the highest level is reached in which interactions between every part of the structure is considered.

The inventive method may be characterized as a full-wave approach. The inventive fill-wave approach is rigorous and therefore more accurate than conventional approximation methods. The computation of each sub-domain is isolated from the others, providing a natural interface to parallel computing. Aided by the power of parallel computers, the novel EPA-based system and method provides for efficient analysis of complex circuit structures. In the case of package modeling, templates are easily created for parameterized analysis. Each subdomain or template is constructed to be reused in a case where repetitive structures such as via arrays or fan-out wiring and very large structures are readily handled, regardless of variations of dimension. 

1. A method for performing a signal integrity analysis on an integrated circuit (IC) comprising a plurality of scatterers by dividing the scatterers into subgroups using a nested Huygens' equivalence principle algorithm, and solving a set of equations realized thereby using a reduced coupling matrix, comprising the steps of: decomposing the IC design into a plurality of small non-overlapping circuit sub-domains, wherein each of the sub-domains is formed as a small, enclosed region; analyzing each sub-domain independent of the other sub-domains using only EFIE to represent the interactions of each sub-domains with the other sub-domains as equivalent currents on equivalent surfaces of the plurality of sub-domains; grouping together neighboring equivalent sub-domains to form larger sub-domains using equivalent currents on equivalent surfaces to represent the interactions of the sub-domains; and repeating the steps of analyzing and grouping the sub-domains until the grouping approaches a box comprising the entire domain, and that the domain interactions between every sub-domain has been analyzed.
 2. The method as set forth in claim 1, wherein the computation of each sub-domain is isolated from the others in order that the method is readily implemented using a parallel computing approach.
 3. The method as set forth in claim 2, implemented as an equivalence principle algorithm (EPA).
 4. The method as set forth in claim 1, wherein the IC design includes an IC package design, and includes step of creating package design templates for parameterized analysis.
 5. The method as set forth in claim 4, wherein each sub-domain or template is reused when comprising a repetitive structure.
 6. The method as set forth in claim 1, wherein simulations may be performed at frequencies up to 10+ GHz.
 7. A systematic mathematical solver for partitioning an electromagnetic design into a set of sub-structures, the set of sub-structures arranged for quicker processing as smaller structures, and solving a set of equations corresponding to the set of sub-structures while maintaining physically correct wave properties of each sub-structure, comprising: generating a closed surface around each sub-structure for analysis; defining the sub-structures such that any interactions between different sub-structures are defined with respect to the sub-structure's closed surfaces; directly mapping an interaction decomposition of the set of sub-structures to a numerical model that defines a set of surface integral equations; and generating a sparse method of moments (MOM) matrix equation from the set of surface integral equations.
 8. The solver as set forth in claim 7, wherein MOM matrix equation is processed to realize a solution to the set of integral equations. 